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ABSTRACT 

We present the results of a simple but robust morphological classification of a statis- 
tically complete sample of 108 of the most X-ray luminous clusters at 0.15 ^ z ^ 0.7 
observed with Chandra. Our aims are to (a) identify the most disturbed massive clus- 
ters to be used as gravitational lenses for studies of the distant universe and as probes 
of particle acceleration mechanisms resulting in non-thermal radio emission, (b) find 
cluster mergers featuring subcluster trajectories that make them suitable for quanti- 
tative analyses of cluster collisions, and (c) constrain the evolution with redshift of 
the cluster merger fraction. Finally, (d) this paper represents the third public release 
of clusters from the MACS sample, adding 24 clusters to the 46 published previously. 
To classify clusters by degree of relaxation, we use the projected offset of the brightest 
cluster galaxy from the peak (or the global centroid) of the X-ray emission as a mea- 
sure of the segregation between the intracluster gas and dark matter, and also perform 
a visual assessment of the optical and X-ray morphology on all scales. Regarding (a), 
we identify ten complex systems likely to have undergone multiple merger events in 
the recent past. Regarding (b), we identify eleven systems likely to be post-collision, 
binary, head-on mergers (BHOMs), as well as another six mergers that are possible 
BHOMs but probably harder to interpret because of non-negligible impact parameters 
and merger axes closer to our line of sight. Regarding (c), we find a highly significant 
increase with redshift in the fraction of morphologically disturbed clusters (and thus 
a clear decrease in the number of fully relaxed systems) starting at z ~ 0.4, in spite of 
a detection bias in our sample against very disturbed systems at high redshift. Since 
our morphological diagnostics are all based on imaging data and thus sensitive to pro- 
jection effects, the measured merger fractions should be considered lower limits and 
our list of mergers incomplete, as we are likely to miss systems forming along axes 
close to our line of sight. A larger sample of clusters with high-quality X-ray data in 
particular at high redshift will be needed to trace the evolutionary history of cluster 
growth and relaxation closer to the primary epoch of cluster formation z ^ 1. 



1 INTRODUCTION 



ray flux and luminosity, and the gas temperature (Randall, 



As the largest gravitationally bound objects in the universe, 
clusters of galaxies represent unique laboratories for studies 
of structure formation and evolution on scales ranging from 
tens to thousands of kpc. In the observationally and theo- 
retically well supported framework of hierarchical structure 
formation clusters grow through a series of mergers. The re- 
sulting negative evolution with lookback time of, e.g., the 
cluster mass function has long been recognized as a power- 



Sarazin Ricker||2002l [Poole et al.|[2007| |Wik et al.||2008 ) 

Unaccounted for, such dramatic deviations from virialization 
are likely to bias the results of cosmological cluster studies. 
This is true in the extreme for work based on measurements 
of the baryon fraction in clusters that rely explicitly on the 
assumptions of spherical symmetry and hydrostatic equilib- 
rium ( [Allen et al.||2004[ [2008 ). 



ful tool to constrain cosmological parameters (e.g., [Eke et al. 



T998)[Borgani et al.|200"T][Vikhlinin et al.|2003l|Mantz et al 



2008 2010). However, the very process underlying this evo- 
lution - cluster merging - is much more complex than the 
physics assumed by many evolutionary models. Speciflcally, 
shocks created in cluster collisions have been demonstrated 
to cause temporary but dramatic increases in the immedi- 
ate observables, i.e., the Comptonization parameter, the X- 



While thus somewhat of a nuisance in the context of 
cosmological work, disturbed clusters are sought-after tar- 
gets for investigations of the dynamics and physics governing 
the interactions between the three main cluster constituents, 
dark matter, diffuse gas, and galaxies. The most extreme 
mergers, head-on collisions of massive clusters of roughly 
equal mass, constitute the most energetic events in the Uni- 
verse, involving the conversion of up to 10^^ erg of potential 
energy into thermal and kinetic energy. As a result, massive 
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mergers provide us with a rare opportunity to study a wide 
range of astrophysical phenomena, such as galaxy evolution 
via ram-pressure stripping, the destruction of cool cores, or 
the acceleration of relativistic particles and the role of mag- 
netic fields in the creation of radio halos and relics. 

In this paper, we attempt to facilitate the selection of 
massive clusters most suitable for cosmological or astrophys- 
ical applications by classifying a statistically complete sam- 
ple of clusters according to relaxation state. The pronounced 
segregation of gaseous and (essentially) coUisionless matter 
during a merger represents a simple and powerful diagnos- 
tic to identify systems that experienced a significant merger 
event fairly recently and are still far from relaxation. Our 
classification scheme uses a combination of optical and X- 
ray diagnostics, and is applied to a sample of the most X-ray 
luminous clusters at redshifts between z — 0.15 and z — 0.7. 
Our goals are threefold: we aim to (a) compile a list of the 
most complex mergers which, although often too challeng- 
ing for dynamical modeling, are likely to be extraordinary 
gravitational telescopes and prime targets for the detection 
of, e.g., extreme gas temperature or giant radio halos and 
relics; (b) identify the most massive, binary head-on cluster 
mergers that feature sufficiently simple trajectories to allow 
quantitative studies of the properties of dark matter and the 
physics of cluster and galaxy evolution; and c) quantify the 
fraction of dynamically disturbed clusters as a function of 
redshift in order to test theoretical and numerical predic- 
tions. 

Throughout this paper we assume the concordance 
ACDM cosmology (Qm 0.3, Ha 0.7, and Hq 70 
km s~^ Mpc~^). 

1.1 Extreme complex mergers 

The most complicated mergers, i.e. systems that have re- 
cently undergone several mergers along different axes, are 
extraordinary observational targets — and can be extraor- 
dinarily difficult to interpret. They are as far from relax- 
ation as possible for a gravitationally bound system, feature 
pronounced substructure in the distribution of dark matter, 
gas, and galaxies, often exhibit dramatic variations in the 
temperature of the intra-cluster medium (ICM), and cannot 
be assumed to be in hydrostatic equilibrium on any scale. 
Clearly these systems grossly violate the assumptions made 
in cosmological models of the growth of structure built solely 
upon spherical gravitational collapse. 

The exceptional properties of active, serial mergers ren- 
der them extremely valuable in other regards though. Ow- 
ing to the large regions of high mass surface density in 
their centres, complex mergers often make exquisite grav- 
itational telescopes whose large Einstein radii allow stud- 
ies of distant background objects over larger areas than ac- 
cessible with relaxed cluster lenses. Examples at intermedi- 



ate redshift include MACS J0717. 5-f3745 ([Ebehng, Barrett 
Donovan] |2004| [Ebeling et al.||2007| [Zitrin et al.||2011|) 



and M ACS J1149.5-h2223 ( [Ebeling et al.||2007| [Smith et al 
|200 9 ) . Although often too complex to be easily modeled 
dynamically, some exceptionally active mergers have been 
used successfully also for astrophysical studies of the clus- 
ter proper, on subjects ranging from ICM heating to galaxy 



been found to feature a very high incidence of radio halos 
and relics, making them excellent targets for studies of par- 
ticle acceleration mechanisms and the role of magnetic fields 
in the confinement of relativistic plasma (see Mayor et al. 
Giovannini et al^[2009[ and references therein). 
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1.2 Binary, head-on cluster mergers (BHOMs) 

Because of their relatively simple geometry, binary head- 
on mergers (henceforth BHOMs) of massive clusters are the 
most promising targets for studies of the physics of mergers 
if a robust physical interpretation of the observational evi- 
dence is sought. Not only does a head-on collision maximize 
the observable signal from all types of interactions between 
the cluster components, the facts that only two systems are 
involved and that a small or negligible impact parameter en- 
sures a linear collision trajectory also greatly reduce system- 
atic uncertainties and facilitate the modeling of the merger 
history of BHOMs. 

In recent studies of cluster mergers, the different dy- 
namical behaviour of the viscous ICM and the essentially 
coUisionless galaxies and dark matter has yielded com- 
pelling observational evidence of the existence of dark mat- 
ter ( Clowe et al.[[2006 Limousin et al.[|2010 Richard et al. 



2010 ) and has, in addition, been used to place upper limits 
on the dark-matter self-interaction cross section, cr/m, of the 
putative dark-matter particle ( jMarkevitch et al.|2004 Clowe 
et al.|2006 i Bradac et al . 200^. BHOMs involving clusters of 
similar mass provide the tightest constraints in this regard, 
and also provide the energetically most extreme environment 
for studies of any interaction between the constituents of the 
merging systems. The absence of complications caused by an 
unknown impact parameter makes head-on massive merg- 
ers particularly suitable also for quantitative studies of the 
collision of the gaseous cluster components which leads to 
shock fronts, collisional heating, and the gradual destruction 



of cool cores ([Markevitch Vikhhni n' 2007 ; Milli on et al. 
2010| jEhlert et al.[[2011[ ). In addition, BHOMs have yielded 
valuable insights into the prevailing mechanisms behind ac- 



celerated galaxy evolution in dense environments (e.g. Ma 
et al.[[2010[ ). 



1.3 Redshift evolution of the merger fraction 

Large, statistically well defined cluster samples are required 
to constrain the evolution of cluster properties over cosmo- 
logical timescales. Of particular interest in this context is the 
rate of cluster formation via mergers, a quantity that affects 
several key statistics, among them the cluster X-ray luminos- 
ity, temperature, and mass functions, and which is a strong 
function of the adopted cosmological world model. Numer- 
ical simulations predict an increase in the merger fraction 
with redshift ([Cohn k White|[2005| [Kay et al.[[2007| [Ettori 



Brighenti|2008 ) , consistent with observations of a decline 
in the number of cool-core clusters with redshift (Vikhlinin 



evolution (e.g., 


Ma et al.|2008 Ma, Ebeling & Barrett |2009 


Ma et al.[|2010 


). Last, but not least, cluster mergers have 



et al.[[2007[ [Santos et aP 2008 ; Samue le et"ar][20lit , as well 
as with expectations from essentially any world model based 
on hierarchical structure formation. 

A low-redshift baseline for such studies is provided 
by work on the Highest X-ray Flux Galaxy Cluster Sam- 



ple (HIFLUGCS, [Reiprich Bohringer[ [2002t , a sample 
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of the X-ray brightest clusters, most of which are nearby 
(<2;> =0.053). Katayama et al. ( 2003| derive a value of 



17.5% for the fraction of dynamically disturbed clusters us- 
ing a subset of 40 clusters. A lower estimate of 12% is found 



in a later study by Sanderson, Edge & Smith ( 2009 ) who use 



all 64 clusters in the HIFLUGCS sample. For X-ray lumi- 
nous clusters at somewhat higher redshifts (2;median=0.23) 
Reiprich & Bohringer as well as Hudson et al. C2010) find 



20% for the same statistic; the comparison of these values is, 
however, complicated by the fact that the two studies differ 
both in their cluster selection criteria and in the diagnostics 
employed to quantify dynamical state (see Section [2]). 

Systematic effects also hamper studies of the evolution 
of the merger fraction at higher redshift, such as the inves- 
tigation of [Samuele et al. (2011) who use the detection of 
strong [OII]A3727 emission in the optical spectra of BCGs 
as a proxy for the presence of a cool core. Although these 
authors report a dramatic decline in the number of such 
line emitters as a function of redshift, their local and dis- 
tant cluster samples differ significantly in their mean X-ray 
luminosity, making a direct comparison difficult. Similarly, 
Maughan et al.^ ( ^2008 ) report an increase at z > 0.5 in the 
number of clusters with large shifts in the centroid of the 
X-ray surface brightness as measured from Chandra obser- 
vations of 115 clusters at 0.1 < z < 1.3, but stress that 
their sample is not statistically complete and their analysis 
subject to "unquantifiable biases". 

A census of merging clusters from a well defined and 
uniformly selected sample spanning a significant redshift 
range that allows a robust determination of the evolution 
of the merger fraction would thus be of considerable inter- 
est. We here attempt to perform this determination as well 
as to provide a list of BHOMs and complex mergers, using a 
statistically nearly complete sample of X-ray luminous clus- 
ters at z = 0.15 — 0.7 and the diagnostics described in the 
following. 



2 DYNAMICAL SUBSTRUCTURE 
DIAGNOSTICS 

Substructure in clusters is notoriously difficult to quantify 
reliably. In principle, spatially resolved X-ray spectroscopy 
can probe the thermodynamics of the ICM in great detail, 
revealing merger-induced shocks and temperature variations 
as well as cold fronts and cool cores, all of which represent 
sensitive diagnostics of a system's relaxation state. In prac- 
tice, however, the required exquisite spatial resolution and 
photon statistics can be attained only for few nearby and 
very X-ray bright clusters. Observationally less "expensive" 
probes of substructure in the ICM use only the projected X- 
ray surface brightness distribution to compute, e.g., power 
ratios (Buote & Tsai 1995) or centroid variations and ax- 



ial ratios ( Mohr et al.|l995| ). More recent studies found the 
combination of data from different wavelength regimes to 
be particularly powerful for the purpose of allowing a crude 
but robust classification. fHudson et al. (2010) investigate 
the usefulness and efficiency of several diagnostics for the 
identification of cool-core clusters, considered dynamically 
fully relaxed systems. Among other criteria, the presence of 
a central temperature drop, short cooling time scales, signif- 
icant mass deposition rate, and the presence of diffuse large- 



scale radio emission are all found to reliably characterise a 
system's dynamical state. However, Hudson and co-workers 
also demonstrate that an even simpler diagnostic, namely 
the projected separation of the locations of the brightest 
cluster galaxy (BCG) and of the peak of the X-ray emission 
from the diffuse intra-cluster gas, constitutes just as robust 



an indicator. Sanderson, Edge & Smith ( 2009 ), in their study 



of the dynamical state of clusters at 0.15<2;<0.3, use two 
different flavours of such a BCG-X-ray offset to identify dy- 
namically disturbed systems: one using the location of the 
X-ray peak, the other using the overall centroid of the diffuse 
X-ray emission. 

The physical reason for which the mentioned X- 
ray /optical offsets allow the straightforward identification of 
disturbed clusters lies in the different dynamical behaviour 
of gas and galaxies during cluster mergers that we mentioned 
before. On cluster scales, galaxies can be taken to be fast- 
moving, pointlike particles that have a very low cross sec- 
tion for collision. By contrast, the intracluster gas is highly 
collisional and, during a head-on cluster merger, will be sub- 
ject to intense ram pressure and shock heating. As a result, 
the galaxy populations of the merging systems will proceed 
nearly unimpeded, slowed only by gravity through the pro- 
cess of dynamical friction, while the strongly self-interacting 
viscous cluster gas is slowed by the collision. Since the re- 
sulting offset will diminish as the combined system gradu- 
ally settles into dynamic and hydrostatic equilibrium, a pro- 
nounced spatial segregation of gas and galaxies constitutes 



powerful evidence of a recent or ongoing merger (e.g.. Smith 
et al.|200"5l Shan et al.|2010 |) ; conversely, excellent alignment 



of the BCG with the intra-cluster gas distribution is typi- 
cal of relaxed systems which often also exhibit significant 



radiative cooling in the core region (Mittal et al. 2009). If 



the merger axis lies approximately in the plane of the sky, 
the apparent separation is maximised, allowing observers to 
deduce the direction of motion and to more readily separate 
the merging subclusters for detailed study. 

How exactly this "X-ray /optical offset" is quantified af- 
fects the respective diagnostic's sensitivity to the complexity 
of the merger (more than two merging components) and/or 
the inclination of the merger axis relative to the plane of 
the sky. Specifically, the observed offset of the BCG from 
the peak of the diffuse X-ray emission is most useful to de- 
rive the projected trajectory of the cores of merging clusters, 
whereas its projected separation from the global centroid of 
the cluster X-ray emission is often more reliable to iden- 
tify morphologically complex systems with significant X-ray 
substructure. We here use both statistics in order to maxi- 
mize our ability to reliably identify disturbed systems and 
mitigate the dependence on merger geometry and viewing 
angle. 

Finally, we use, as a third, qualitative diagnostic of the 
relaxation state of all clusters the result of a visual classi- 
fication of the morphology of the X-ray contours and the 
degree of alignment between the BCG and the X-ray peak, 
using the simple classification scheme devised by |Ebeling| 



et al.| ( |2007p . The assigned codes and classification criteria 



(from relaxed to extremely disturbed) are 1 (pronounced 
cool core, perfect alignment of X-ray peak and single BCG), 
2 (good X-ray/optical alignment, concentric contours), 3 
(non-concentric contours, obvious small-scale substructure), 
and 4 (poor X-ray/optical alignment, multiple peaks, am- 
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biguous BCG). We estimate the uncertainty in the assigned 
morphological code to be less than 1. 



3 CLUSTER SAMPLE 

Our sample is selected from a master list of clusters iden- 
tified in the course of the Massive Cluster Survey (MACS; 
Ebeling, Edge k Henry] [200T| [Ebeling et aL][2Q07l [Mof. 
While MACS, by definition, is limited to the most distant 
systems (z ^ 0.3), clusters (or candidate clusters) at lower 
redshift were also flagged by the MACS team during their 
screening of thousands of X-ray sources from the Bright 
Source Catalogue compiled in the course of the ROSAT All- 
Sky Survey (RASS, |Truemper||1993| . We here draw from 
this master list of RASS clusters with /rass ^ 1 x 10~^^ 
erg s""*^ cm~^ (0.1-2.4 keV), a database that covers the en- 
tire extragalactic sky observable from Mauna Kea, i.e., the 
solid angle within the boundaries defined by \b\ > 20° and 
-40° < ^(J2000) < 80°. 

Application of the additional selection criteria z > 0.15 
and Lx,RASS >5 x lO'^^ erg s"^ (0.1-2.4 keV) yields a list 
of 129 clusters (Fig.[T]). The quoted lower redshift limit was 
imposed to ensure that the diffuse X-ray emission from all 
clusters is well contained within the field of view of the ACIS 
(Advanced CCD Imaging Spectrometer) instrument aboard 
the Chandra X-Ray Observatory (CXO). Carrying the only 
X-ray telescope that provides sub-arcsec resolution, Chan- 
dra is the observatory of choice for the unambiguous dis- 
tinction between small-scale structure (especially cold cores) 
and X-ray point sources, an essential requirement for a ro- 
bust assessment of the X-ray morphology of the clusters un- 
der stud}|^ The quoted luminosity cut was applied to select 
an unbiased sample of similarly X-ray luminous (and thus, 
by inference, massive) clusters across the full redshift range, 
and to increase the probability of the selected systems hav- 
ing been observed by CXO. As shown in Fig. [l] our cluster 
sample is volume complete out to z ^ 0.4 and X-ray flux 
limited beyond this redshift. Of the 129 clusters, 108 have 
been observed with ACIS as of July 2011. 



4 DATA AND ANALYSIS 
4.1 X-ray Data 

Table 1 lists the exposure time for each cluster observation. 
We merged data from multiple observations wherever doing 
so increased the combined exposure time by at least 30%. 
We used CI AO 3.3.0 to create images in the 0.5-7 keV band, 
which we then smoothed using the adaptive- kernel algorithm 



ASMOOTH ( [Ebeling, White k Rangarajan||2006 ). 

The only quantitative X-ray characteristics required for 
our study are the locations of the peak and of the global 
centroid of the emission. We calculate the former from the 
adaptively smoothed image to mitigate the impact of fluctu- 
ations in the photon statistics, in particular for observations 



of short duration. In order to measure positions to better 
than 5kpc accuracjj^we compute the peak position as the 
photon-weighted mean within a circle of 2'' radius around 
the brightest X-ray pixel. We confirmed that the result of 
this process is robust and insensitive to the precise choice of 
radius. We also stress that, thanks to the sub-arcsec angu- 
lar resolution of CXO, X-ray point sources at or very near 
the location of the peak of the diffuse X-ray emission can 
be masked out efficiently and have no impact on the mea- 
surement of the position of the peak of the diffuse emission. 
Where multiple peaks of comparable brightness are encoun- 
tered (peak brightness difference of less than a factor of two) 
we adopt the brighter peak for the computation of the X- 
ray/BCG offset but flag the respective cluster for separate 
visual assessment. 

The second X-ray characteristic required for this work 
is the location of the overall centroid of the diffuse X-ray 
emission. We compute this position as the photon- weighted 
mean of the coordinates of all pixels exceeding a surface- 
brightness threshold of 1.2 x 10~^ counts s~^ for ACIS-I 
data, and 1.6 x 10~^ counts s~^ for ACIS-S, thereby account- 
ing for the higher instrumental background of the back- 
illuminated CCD at the ACIS-S aim point. The quoted 
surface-brightness thresholds also mark the count rate lev- 
els of the lowest contours in the X-ray/optical overlays used 
to assess the overall cluster morphology (see Fig. |2] for an 
example) . 

Our ability to discern multiple X-ray peaks (or, for that 
matter, to accurately measure the position of any X-ray 
peak) depends on the depth of the respective X-ray image, 
as well as on the X-ray brightness of the cluster and its X- 
ray morphology, all of which vary significantly between clus- 
ters. Figure [2] illustrates this point by showing, side by side, 
the adaptively smoothed images of A 2744 (MACS j0014.3- 
3022) from the 74 ks observation obtained with ACIS-I, as 
well as from just the first 1.5 ks of the same data set. As is ap- 
parent from the comparison, the global centroid of the X-ray 
emission is hardly affected by the depth of the observation 
used; the exact location of the X-ray peak(s), however, can 
change significantly as a function of photon statistics. To 
assess whether variations in the photon statistics within our 
sample affect our ability to obtain a robust morphological 
classification at all redshifts we show in Fig. |3] the photon 
count (net of background) for all clusters as a function of 
redshift. As a result of observer bia^ photon statistics are 
poorest at roughly the midpoint of the redshift range of our 
sample. Although still sufficient for reliable determinations 
of the X-ray centroids, the CXO photon statistics for some 
clusters at 2:^0.45 may, for certain morphologies and view- 
ing angles, be too poor to allow accurate measurements of 
the X-ray peak location. 



^ The original detections (median, net photon count of 54) from 
the RASS allow only a determination of the total flux / luminos- 
ity, and even for the brightest system (574 net photons within 
the detect cell) even a crude morphological classification is made 
impossible by the poor angular resolution of the RASS. 



^ The pixel scale of ACIS, 0.5^^ corresponds to 6 kpc at z = 0.68, 
the highest redshift of all clusters in our sample. 
^ The most distant X-ray luminous systems receive special atten- 
tion (and thus more CXO exposure time) owing to their impor- 
tance for cosmological studies. 
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4.2 Optical Data 

In addition to X-ray images of good angular resolution (for 
the computation of the X-ray peak and centroid positions), 
we need optical images that allow the unambiguous identi- 
fication of the brightest cluster galaxy (BCG). High-quality 
colour images of all clusters at z > 0.3 in our sample were 
obtained with the UH2.2m telescope in the course of the 
compilation of the MACS sample. At lower redshift, we per- 
formed similar observations for a subset of our sample and 



rely on SDSS (Sloan Digital Sky Survey, Stoughton et al. 



2002 ) and DSS (Digitized Sky Survey) data for the remain- 



der. A brief description of each of these imaging datasets 
follows. 

UH2.2m Data: Imaging with the UH2.2m telescope was 
performed for 77 of the 108 clusters in our sample between 
1999 and 2010, using the Tek2048 CCD (pixel size 0.22'') 
and 240 second dithered integrations in the V ^ R and / 
filters. Obtained in photometric conditions and sub-arcsec 
seeing, these images allow the determination of the BCG 
centre with negligibly small errors. 

SDSS Data: For 23 clusters SDSS colour images were 
produced from data taken through the ^, r, and i filters 
(Fukugita et al.(1996). Although not as deep as the UH2.2m 



images and also featuring worse angular resolution (median 
seeing ~ 1.4'' and using 0.5" pixels; Abazajian et al.|2009 ), 
the SDSS images still yield BCG positions with negligibly 
small errors. 

HST Data: Three clusters (A 2146, A 2163, and A 2667) 
located outside the SDSS footprint and without UH2.2m 
images have been observed with the Hubble Space Telescope 
(HST) allowing highly accurate measurements of the BCG 
position. 

DSS Data: For the remaining five clusters, the only 
readily available optical data hail from the DSS, which is 
much shallowei^ of poorer angular resolutior[^ and of in- 
ferior cosmetic qualitjj^ than the aforementioned imaging 
databases. We use infrared, red, and blue images from the 
DSS to create colour images. The relatively poor quality of 
the DSS data leads to significant uncertainties in the deter- 
mination of the BCG centre. 



4.3 Astrometric Corrections 

The primary diagnostics used by us to characterise the dy- 
namical state of clusters and to identify merging systems 
are the offsets between the location of the BCG and that 
of the X-ray peak or global X-ray centroid, respectively. Ei- 
ther diagnostic thus mandates that our optical and X-ray 
images be astrometrically aligned as accurately as possible. 
The magnitude of the required astrometric corrections can 
be estimated from the absolute astrometric accuracy of the 



For instance, the limiting magnitude of the DSS is ~ 21.5 
(|Gal et al.||2004| > compared to ^ ~ 23 for the SDSS ( [York et al7| 
2000|. 



^ I'he median seeing is ~ l.h" and the pixel size 1" . 
^ Much of the DSS data was obtained from scans of photographic 
plates, a process that, in addition to objects of astrophysical in- 
terest, also faithfully recorded fingerprints, scratches, and other 
artefacts. 



datasets involved. On the X-ray side, CXO provides astro- 
metric solutions that are accurate to better than 0.8" within 
3 arcmin of the ACIS aimpoinl[^ For our optical data, we can 
expect astrometric uncertainties for our UH2.2m images of 
less than 0.5" after processin g with WCSToo ls ( |Mink|2006| ), 
less than 0.2" for the SDSS ( |Pier et al.|2003|, and less than 



0.4" for the DSS images used (Lasker et al 



2008). The re- 



sulting error in the relative astrometry can thus reach, and 
occasionally even exceed, l", which corresponds to 5 kpc at 
the median redshift [z — 0.35) of our sample - a significant 
contribution to the overall error budget. 

In order to eliminate, or at least minimise, this source 
of measurement error, we compute corrections in the relative 
astrometry between the optical and X-ray images by identi- 
fying the optical counterparts to the brightest X-ray point 
sources within the field of view common to both images. 
Typically, over half a dozen matches with X-ray/optical off- 
sets that are consistent to within 0.3" are found for each 
cluster. Figure^ shows a histogram of the astrometric cor- 
rections applied; the distribution is consistent with the ex- 
pectations discussed above. 



4.4 Merger classification criteria 

While the morphological diagnostics discussed in Section |2] 
are well suited to allow us to select the most disturbed sys- 
tems, several selection biases enter. 

Most importantly, morphological classification based 
solely on imaging data is subject to projection effects. Given 
that all clusters selected for this study are intrinsically X-ray 
luminous, projections involving physically unrelated fore- or 
background structure are improbable. Hence, overestimates 
of the complexity of a system, and thus of the merger frac- 
tion, are negligible. The opposite effect, however, is bound 
to be present, as mergers proceeding along or close to our 
line of sight are likely to go unrecognized. We shall return 
to this important point when discussing the results of our 
study. An additional complication arises from the fact that 
we are interested only in mergers observed after the primary 
collision, thereby rendering them useful laboratories for in- 
depth study of the interactions of all three cluster compo- 
nents (Sections |1.1[ |1.2| . Visual scrutiny of the combined 
X-ray/optical appearance of a cluster is usually sufficient 
to eliminate close double clusters or pre-collision mergers, 
but - alas - again not for prolate systems forming along 
an axis close to our line of sight. Finally, the subdivision of 
disturbed clusters into complex mergers and BHOMs is not 
always possible and necessarily subjective to some extent. 

These challenges notwithstanding, we attempt to iden- 
tify complex mergers as well as BHOM candidates, all of 
them observed after the primary collision, by using the fol- 
lowing criteria: 

Morphological class: Any ongoing or recent merger 
along an axis sufficiently misaligned with our line of sight 
will have a value of 3 or 4 for the morphological code de- 
scribed in Section[2] A high value for the morphological class 
is thus a necessary but not a sufficient criterion for the se- 
lection of both complex mergers and post-collision BHOMs; 



See http://cxc.harvard.edu/cal/ASPECT/celmon/ 
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X-ray-centroid/BCG separation: A large offset be- 
tween the global centroid of the X-ray emission and the lo- 
cation of the BCG may be an indicator of a post-collision 
merger in which the ICM of the merging clusters have coa- 
lesced into a single central gas cloud; however, large offsets 
are also observed for pre-collision mergers. Conversely, view- 
ing angle and merger phase may conspire to yield a small 
X-ray-centroid/BCG separation even for very disturbed sys- 
tems. A high value for this diagnostic is thus (with the noted 
caveat concerning pre-collision mergers) a sufficient but not 
a necessary selection criterion; 

X-ray-peak/BCG separation: A large separation be- 
tween the X-ray peak and the location of the BCG is an 
unambiguous sign of an ongoing merger. Correctly inter- 
preting a small separation can, however, be a difficult task 
since this metric depends on both the stage of the merger 
process (BCG and surviving cool core will align again ap- 
proximately at turnaround) and the impact parameter (in 
particular for mergers of clusters of significantly different 
mass). A high value for the X-ray-peak/BCG separation is 
thus again a sufficient but not a necessary selection criterion 
for our purposes. 

In view of the complexity and enormous variety in the 
projected geometry and dynamics of real-life cluster merg- 
ers, we therefore require likely post-collision mergers to a) 
meet our morphological selection criterion (class 3 or 4) and 
b) feature a large offset of the BCG from either the X-ray 
centroid or the X-ray peak. 

Noticeably absent from the discussion above is a cri- 
terion that clearly distinguishes between complex mergers 
and BHOMs, i.e., a criterion that assesses whether two, and 
only two, clusters are involved and clearly discernible in the 
optical image of a BHOM candidate. Ideally we would hence 
like to add a fourth criterion akin to: 

BCG-BCG separation: Unless it proceeds along an axis 
very close to our line of sight, a two-cluster merger observed 
at any stage other than core passage should feature two 
well separated BCGs, one for each of the two participating 
clusters; for a post-collision BHOM, the X-ray centroid will 
be located between them. 

Unfortunately, we found this BHOM selection criterion (in- 
tuitive and straightforward as it seems) to be almost im- 
possible to implement in practice, the primary reason being 
that, for mergers of very different mass, the two brightest 
galaxies within the lowest X-ray contour are frequently both 
members of the more massive merger component. This is 
true, for instance, for the Bullet Cluster: the two brightest 
galaxies are both members of the (unrelaxed) main cluster. 
While this criterion can thus not be applied without visual 
scrutiny of the respective system, it is useful in principle, as 
illustrated by the case of the highly complex merger A 2744 
which features no fewer than five galaxies of essentially equal 
brightness (within the measurement uncertainties. Fig. |2|. 

In recognition of the difficulties inherent in the subclas- 
sification of merger types we therefore apply the first three 
criteria to select extreme mergers of all kinds, and then re- 
sort to visual inspection of all candidates in order to dis- 



5 RESULTS 

5.1 X-ray /optical separation and morphological 
code 

Table 1 lists the values of all three quantitative diagnos- 
tics used in this study, i.e., the assessed global morphology 
described in Section |2] as well as the projected metric sepa- 
rations between the BCG and the X-ray peak and centroid, 
respectively, for each of our 108 clusters. Also listed are the 
la errors of the latter two measurements. We find the distri- 
bution to be roughly log-normal for either type of separation 
(Fig. [5|, and to be centred at 11.5 and 21.2 kpc for the off- 
set of the BCG from the X-ray peak and from the X-ray 
centroid, respectively. 

15% (16 out of 108) of the clusters in our sample have 
projected X-ray-peak/BCG separations of at least 50 kpc, 
slightly higher than, but (within the uncertainties) consis- 
tent with, the percentage of 12% (8 out of 64) reported by 
[Hudson et al.| ( |2010|) for the HIFLUGCS sample. Application 
of the same cut (>50 kpc) to the X-ray-centroid/BCG sep- 
arations selects 30% (32 out of 108) compared to a value of 
20% (13 out of 65) found for the LoCuSS sample (Sander- 



son, Edge & Smith 2009). Both of these comparisons are. 



however, problematic because of the fact that HIFLUGCS 
is X-ray flux limited, and LoCuSS uses Galactic nu as a 
selection criterion, whereas the sample used in this work is 
primarily X-ray luminosity limited. In addition, the redshift 
ranges covered by the three cluster samples differ substan- 
tially (HIFLUGCS: <z> = 0.05; LoCuSS: <z> = 0.23; this 
study: <z> = 0.35). 



tinguish between complex mergers (Section 5.2.1 ) and likely 
post-collision BHOMs (Section [5.2.21). 



5.2 Extreme and active mergers 

A total of 29 clusters in our sample meet the criteria dis- 
cussed in Section [4.4| i.e., they have been assigned a mor- 
phological code of 3 or 4 and feature a large offset of the 
BCG from the X-ray centroid or the X-ray peak. We imple- 
ment the latter criterion by requiring the respective offset 
to exceed a threshold value of 71 or 42 kpc, i.e., the mean of 
the respective distribution of separations plus one standard 
deviation, rounded to 1 kpc precision (Fig.[5]). 

All systems thus selected are far from relaxed and thus 
prime candidates for studies of the physics of cluster merg- 
ers; however, not all of them are BHOMs, and some may 
not meet the important requirement that we observe them 
after the primary core passage. As discussed in Section [4. 4[ 
it is difficult to devise an objective, quantitative criterion 
for either char act istic. We therefore screen all candidates vi- 
sually, focusing on the systems' appearance in the optical 
and X-ray waveband. Our assessment then uses primarily 
the merger's morphological complexity as well as the direc- 
tion of motion of the merger participants inferred from the 
BCG-X-ray-peak offsets. 

As a first result of this screening process we elimi- 
nate MACS J0159.0-3412 and ZwCl 1459.4+4240 because, 
in both cases, the cluster emission falls so close to the ACIS-I 
chip gaps as to raise doubts about the reliability of the X-ray 
morphological classification listed in Table 1. 
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5.2.1 The most complex massive mergers 

Scrutiny of the optical and X-ray appearance of the 27 ini- 
tial candidates leads to the classification of 10 systems as 
complex mergers that have recently undergone, or are still 
undergoing, multiple mergers along differing axes. We briefly 
discuss them individually below. 

A 2744^ a well studied cluster, features multiple galaxies 
of very similar brightness, rendering the identification of 
a BCG difficult. Although A 2744 could be mistaken for a 
BHOM when only its X-ray appearance is considered, it fea- 
tures X-ray/optical offsets that are inconsistent with the lin- 
ear trajectories expected in such a scenario. Indeed, detailed 
analyses at optical and X-ray wavelengths showed the sys- 
tem to be the result of a recent line-of-sight primary merger 
combined with a secondary merger proceeding at high im- 
pact parameter ( Kempner David|2004 Owers et al.|2011 



[Me rten et al. 2011). Consistent with such an active dy- 
namical history, A 2744 is also known to host a radio halo 
( Boschin et"aL]|2006D . 

A 520 is a complex system whose unusual morphology led to 
claims of the presence of a "dark core" (Mahdavi et al. 2007|). 
A dynamical analysis of the galaxies in the field suggests 
a multiple merger at the intersection of three large-scale 
filaments ( Girardi et al.|200'8 ). Evidence of past and ongoing 
merging is found in the form of a bow shock (one of less than 
a handful of supersonic shocks detected so far in cluster 
mergers) and a complex radio halo ( [Markevitch et al.|2005 ). 
MACS J0717. 5+3745 is the most disturbed massive cluster 



known at z>0.5 (Ebeling et al. 2007), accreting matter along 
a 6 Mpc long filament (Ebeling, Barrett k Donov an|2004| . A 
detailed study of the ICM and the cluster galaxy population 
yielded compelling evidence of a triple merger, resulting in 
ICM temperatures as high as 20 keV and galaxy transforma- 
tions caused by ram pressure ( |Ma et Ebeling 
k Barrett||2009| |Ma Eb eling 2011). Shocks triggered by 



a series of mergers are likely responsible for the accelera- 
tion of electrons in the ICM to relativistic energies, making 
MACS J071 7. 5+3745 a source of extreme non-thermal emis- 
sion at radio wavelengths (Edge et al.|2003 [Bonafede et al^ 



2009 van Weeren et al.|2009| . The system was also success- 



fully targeted in Sunyaev-Zel'dovich observations (LaRoque 



et al. 2006). 



A 665 lacks optical evidence of multiple cluster cores in spite 
of a very disturbed morphology. A joint X-ray /optical anal- 
ysis suggests a head-on merger along the line of sight, ob- 
served close to core passage (Gomez, Hughes & Birkinshaw 
200Qt , a scenario that would explain the system's morphol- 



ogy as viewed in projection as well as the observed unimodal, 
but very broad (cr= 1400 km s~^) distribution of radial ve- 
locities for the cluster galaxy population. A 665 is thus an 



example of a likely BHOM (see Section 5.2.2) missed by the 
set of diagnostics used by us here because the merger axis is 
nearly aligned with our line of sight. This merger too is giv- 
ing rise to non-thermal emission in the form of a well studied 
synchrotron radio halo ( Markevitch Vikhlinin]|2001 Fer- 



etti et al.|[2004||Vacca et al.|2010| )" 

A 1300 is a highly disturbed cluster whose optical complexity 
rivals that of A 2744 in that it features no fewer than five 
galaxies vying for the title of BCG. The cluster's complex 
X-ray morphology suggests a series of recent mergers, as well 
as an ongoing merger refiected by the potential presence of a 



cold front. Hosting both a giant radio halo and a radio relic, 
this system features frequently in studies of the connection 
between cluster mergers and diffuse radio emission ( Reid| 



et al.||1999| [Cassano et"ar]|2008[ [2010) but has, to the best 
of our knowledge, yet to attract similar attention at other 
wavelengths. 

MACS JII49. 5+2223 is part of the complete subsample of 
MACS clusters at z > 0.5. Deep follow-up observations 
with HST/ACS revealed spectacular strong- lensing features 
which allowed detailed modelling of the cluster mass dis- 
tribution. Comprising at least four large-scale dark-matter 
halos, MACS Jl 149.5+2223 is one of the most complex clus- 
ters presently known ( Smith et al.|[2009) . The system's ex- 
treme velocity dispersion of 1800 km s~^ (Ebeling et al 



2007) strongly suggests merger activity along our line of 
sight. 

MACS J2129. 4-0741 is the third cluster from the z > 0.5 
subset of the MACS sample to make our list of complex 
mergers and, much like A 665, does so only by virtue of a 
significant offset of the BCG from the large-scale centroid of 
the X-ray emission recorded by Chandra. Featuring the same 
velocity dispersion of 1400 km s"\ MACS J2129.4-0741 may 
be dynamically similar to A 665 also in terms of its recent 
merger history. 

RXJ2228. 5+2036 was discove red during t he ROSAT 
Brightest Cluster Survey (Ebeling et al.||1998| BCS,) and 



followed up in pointed X-ray observations with the ROSAT 
HRI and XMM-Newton. A high mass in excess of 10^^ 
M0 was estimated from X-ray and Sunyaev-Zel'dovich effect 
observations under the assumptions of sphericity and hy- 
drostatic equilibrium ( Pointecouteau et al.||2002l Jia et al. 



2008| ), both of which are unlikely to apply given the clus- 



ter's X-ray/optical appearance (Fig. 10 see also Maughan 
20081 . Although assification of this system as 



et al. 



a merger is seemingly called into question by the absence 
of diffuse radio emission (Venturi et al. 2008), the strong 



correlation between mergers and non-thermal radio emis- 
sion is not without exceptions, as demonstrated by |Russell| 
et al.|(|2011|) whose GMRT observations of A 2146 (see Sec- 



tion 



5.2.2) also failed to detect any sign of diffuse emission. 



MACS J2243. 3-0935 is a highly disturbed cluster according 
to all three of our diagnostics. The merger axis suggested 
by its highly elongated X-ray emission is significantly mis- 
aligned with the axis implied by the two main galaxy con- 
centrations, thus all but ruling out a BHOM scenari o. The 
system's high velocity dispersion of over 1500 km s"^ dEbel- 
ing et ar||2010 ) may be refiective not only of this cluster's 
high mass but also of line-of-sight merger events in its recent 
dynamical history. 

A 2631 exhibits an X-ray /optical morphology similar to 
MACS J2243. 3-0935 but at much lower galaxy surface den- 
sity. Although clearly highly unrelaxed, the cluster is listed 
as a non-detection by Venturi et al. ( 2008 ) who conducted 
a search for diffuse radio emission in 50 X-ray selected clus- 
ters. The absence of non-thermal radio emission from this 
and several other obvious mergers presents a formidable 
challenge for theories of particle acceleration during mas- 
sive cluster mergers. 
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5.2.2 BHOMs 

The remaining 17 clusters are all potential or probable early 
post-collision BHOMs, and are discussed individually below 
(in order of right ascension). 

Primary candidates 

MACS J0025. 4-1222 is a textbook case of the optical and 
X-ray morphology expected for a post-collision BHOM. It is 
also one of only two systems studied already in depth in at- 
tempts to constrain the collisional properties of dark matter 
(the other one being the Bullet Cluster). The simple geome- 
try of the merger and the wide separation of its components 
made MACS J0025.4-1222 a rewarding target also for stud- 
ies of the impact of mergers on galaxy evolution. We refer to 
the literature for a detailed discussion and results ( (Ebeling 
et al.||2007||Bradac et al.|2008|[Ma et al.|2Q10t . 
MACS JOI40. 0-0555 is very similar in appearance to 
MACS J0025.4-1222, suggesting a post-collision BHOM in- 
volving two clusters of comparable mass. Although the met- 
ric separation of the two components of about 250 kpc may 
pose challenges for the separation of the corresponding dark- 
matter halos in a weak-lensing study, the likely high mass of 
either cluster makes this BHOM an excellent candidate for 
constraining dark-matter properties following the approach 



taken for MACS J0025.4-1222 (Bradac et al. 2008). A first 



characterization of the system based on all existing optical 
and X-ray data is given by Ho et al. (2011, in preparation). 
MACS JO417. 5-1154 has a comet-like appearance in X-rays, 
with an extremely X-ray bright and compact core that is co- 
incident with the optical core of one of the two merging clus- 
ters, trailed by diffuse emission towards, and encompassing, 
the second cluster. The system's X-ray morphology strongly 
suggests a high-velocity encounter akin to the scenario ob- 
served in A 2146 or the Bullet Cluster. The absence of a 
noticeable offset between the locations of the BCG and X- 
ray peak of the compact core suggest that the merger is 
observed either close to turnaround and / or proceeds along 
a merger axis that is significantly tilted in the direction of 
our line of sight. A detailed analysis of the apparent shock 



front to the SW of the compact core (cf. Fig. 11 ), as well as 



of the mass distribution of both components as revealed by 
strong- and weak-lensing observations, is underway. 
A 521 deserves special mention: although the main cluster 
in this merger appears morphologically unremarkabl^ its 
extremely low surface density, both in X-rays and in terms 
of projected galaxy numbers, make it stand out among clus- 
ters of similar X-ray luminosity at this redshift (z=0.253). 



The presence of an extended radio source ( Giacintucci et al. 



2006 Dallacasa et al. 2009) as well as of multiple compo- 



nents in the redshift distribution of cluster galaxies (Ferrari 
et al.|2003 ) strongly suggests a long history of merger events, 
possibly along different axes. In spite of the complexity of 
the primary cluster component, the ongoing merger that ap- 
pears to occur close to the plane of the sky makes A 521 a 
promising target for quantitative studies of the physics of 
cluster mergers ( [Ferrari et al.||2006 ). 



^ Note that the X-ray appearance of this system in Fig.[lljis mis- 
leading due to the ACIS-I chip gaps truncating the very extended 
emission from the presumably more massive southern merger 
component. 



MACS J0553. 4-3342 is again a BHOM that is clearly ob- 
served past core passage and that appears to involve at least 
two similarly massive clusters. The large metric separation 
of the involved clusters makes it likely that the merger axis 
lies approximately in the plane of the sky. Whether a third, 
less massive cluster component is involved (as suggested by 
both the ICM and galaxy distributions), and whether the 
merger is viewed sufficiently shortly after the initial col- 
lision to allow the participating dark-matter halos to still 
be discernibly separated will be revealed by an ongoing X- 
ray/lensing study using deep Chandra and HST/ACS ob- 
servations. 

MACS J1006. 9+3200 resembles A 1682 (see below) in its X- 
ray /optical morphology. Clearly observed after the primary 
collision (possibly one of several recent merger events), the 
system may be near turnaround, as suggested by the wide 
separation of the two subclusters and the direction of the 
BCG/X-ray peak offset of the minority component. To the 
best of our knowledge, no in-depth study of this merger has 
been conducted to date. 

A 1682 appears at first glance like a likely pre-collision 
merger of two systems of very different mass (Fig. 11). A 



pair of BCGs in the primary cluster and diffuse X-ray emis- 
sion between the merging systems, however, constitutes con- 
vincing evidence of prior interaction. This interpretation is 
supported by the detection of very complex diffuse radio 
emission, interpreted as a halo and two relics by |Venturi| 



et al. 



(2008). 



A 1 758 is an obvious post-collision BHOM (prominent X-ray 
emission between the two participating clusters) in which 
both merger components appear to have retained some of 
their original cold cores, visible as compact X-ray peaks 
near the respective BCG. An existing weak-lensing study 
of A 1758 based on ground based optical data (Okab e IE] 
Umetsu|2008 ) finds two well separated mass concentrations 



and assigns the system a total mass of (5± 1.5) X 10 M(: 
and a velocity dispersion of (670 ± 110) km s~^. Although 
thus likely less massive than MACS J0025.4-1222, A 1758 
may be another good candidate for a measurement of a/m. 
The system is of interest also as a future merger on much 
larger scales, given the presence of a second cluster, of lower 
X-ray luminosity but at the same redshift, about 2 Mpc 
south of the system discussed here ( Rizza et al.|1998 David 
& Kempner 2004). Observations conducted with the Spitzer 
observatory revealed an unusually high rate of star forma- 
tion, primarily from the northern, main cluster, with much 
of the presumably merger-induced activity originating in ob- 
scured galaxies near the cluster core and at the extreme out- 
skirts (Haines et al.|2009|. 

A 1914 is a well known merger with a morphologically very 
complex core, but qualifies for "post-collision BHOM" sta- 
tus by virtue of the presence of two distinct (if close) optical 
cluster cores on either side of both the peak and the over- 
all centroid of the X-ray emission. However, the noticeable 
misalignment between the elongated X-ray contours in the 
core and the merger axis suggested by the line connecting 
the BCGs of the two components raises the possibility that 
the collision occurred with a significant impact parameter, 
rather than head-on. In addition, the large extent and near- 
sphericity of the X-ray emission on larger scales may be 
the result of the merger axis being distinctly inclined with 
respect to the plane of the sky. This interpretation is sup- 
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ported by a large discrepancy between mass estimates based 
on X-ray and weak-lensing analyses ( Okabe Umetsu|2008 



Zhang et al.||2010 ). Featuring a powerful radio halo, A 1914 
is also part of a larger cluster sample for which G ovoni et al.| 
(2004) investigated the relation between variations in the 



ICM temperature and the presence and location of diffuse 
radio emission. 

A 2146 is almost a clone of MACS J04 1 7.5-1154 in terms 
of optical and X-ray morphology. A recent in-depth study 



of a deep CXO/ACIS-S observation of this system (Russell 
et al.|[20To| finds two pronounced shock fronts (with Mach 



numbers of 2.2 and 1.7, respectively), rendering this merger 
the closest rival yet to the Bullet Cluster. Although the es- 
timated velocity of the bullet in A 2146 is lower than the 
one measured for IE 0657-56, the ongoing analysis of this 
BHOM is likely to yield competitive constraints on a/m. 
The system is also a prime candidate for a search for evi- 
dence of accelerated galaxy evolution as a result of the col- 
lision. Surprisingly for such an extreme merger, and so far 
inexplicably, A 2146 shows no signs of diffuse radio emission 
( Russell et al]|2011 ). 

MACS J1731. 6+2252 exhibits an X-ray /optical morphology 
similar to that of A 2744 or the Bullet Cluster (although at 
low fidelity, owing to the short exposure time of the Chandra 
observation), in that a merger of two clusters of very differ- 
ent mass is viewed well after collision. The main component 
of MACS J1731. 6+2252 shows no significant X-ray/optical 
offset though, suggesting either a collision at large impact 
parameter and/or a significant inclination of the merger axis 
with respect to the plane of the sky. Like A 2744, this system 
may also be a merger of more than two components. A much 
longer Chandra observation would allow a better estimate 
of the three-dimensional merger trajectories to be obtained. 

Secondary candidates 

The following systems show compelling evidence of be- 
ing post-collision BHOMs but feature (projected) separa- 
tions between the BCGs of the merging clusters of less than 
200 kpc, either because the respective mergers proceed along 
axes close to our line of sight, or because the merging com- 
ponents are not yet (or no longer) sufficiently far apart. If 
the apparent close separation of the latter is indeed due to 
a greatly inclined merger axis, projection effects and the 
degeneracy between peculiar velocity and Hubble fiow are 
likely to impede if not prevent the determination of the true 
three-dimensional merger geometry and the physical inter- 
pretation of the observables. If, on the other hand, these 
mergers proceed approximately in the plane of the sky and 
simply happen to be observed by us at a time of small sep- 
aration of the components, all of the systems listed below 
will still be of great interest for studies of all merger-induced 
physical effects, with the likely exception of measurements 
of a/m. 

A 2813 appears moderately relaxed on large scales but ex- 
hibits a double core in both optical and X-ray images Its 
morphology is most easily explained by a merger proceed- 
ing along an axis that is highly inclined with respect to the 



plane of the sky; if so, it is conceivable that we observe the 
system before the primary collision. 

MACS J0358. 8-2955 was only recently recognized as a very 
X-ray luminous cluster z — 0.434 (Ebeling et al. 2010) 
after |Bohringer et al. (20041 misidentified the system as 
A 3192 at z = 0.169. Both at X-ray and optical wavelengths 
two cluster components are clearly separated, but as evi- 
dent are signs of previous interactions. The combination of 
a small projected separation of the cluster cores and an X- 
ray /optical offset of the more compact component that is 
not aligned with the apparent merger axis suggests a col- 
lision at significant impact parameter and possibly a line- 
of-sight merger. A detailed study with Chandra and HST 
is underway in an attempt to better constrain the merger 
geometry. 

MACS J0404-6+1109 features a very interesting morphol- 
ogy with an X-ray peak located between two almost equally 
bright BCGs. The very small separation of the two BCGs 
of only about 100 kpc suggests either a line-of-sight merger 
or a very recent collision. A deeper Chandra observation is 
needed to improve upon the presently poor photon statistics 
and allow a better assessment of the relative distributions 
of coUisional and collisionless matter in the cluster core. 

MACS JO4I6.I-24O3 would easily make our list of primary 
BHOM candidates if the separation of the system's BCGs 
did not fall just shy of the 200 kpc threshold. Illustrating 
an element of subjectivity of our classification, this merger 
(clearly observed after the primary collision, and probably 
even after turnaround) should make an excellent target for 
studies of the physics of cluster mergers. 

ZwCl 0847.2+3617, aka RX J0850.1+3604, is an X-ray dis- 



covery from the Brightest Cluster Sample (Ebeling et al 



1998) that exhibits the classical signature of a post-collision 
BHOM: a single X-ray peak between two galaxy concentra- 
tions. Its other morphological characteristics, however, raise 
doubts about the system's usefulness for quantitative phys- 
ical investigations. Specifically, the elongation of its X-ray 
emission on large scale runs at right angles to the presumed 
merger axis suggests a complex merger history. In addition, 
the small separation of the mentioned two galaxy concentra- 
tion relative to the extent of the X-ray emission would be 
most easily explained by a merger axis close to the line of 
sight. Diffuse radio emission, as expected from such a mas- 



sive merger, was detected by Rudnick & Lemmerman (2009 ) 
who, however, attribute the emission to a radio galaxy. 

MACS J1115. 2+5320 resembles MACS J0140.0-0555 (and 
thus also MACS J0025.4-1222) but features a smaller sepa- 
ration of the two BCGs. While a partial alignment with our 
line of sight is likely to contribute to the small separation 
observed in projection |Ebeling et al. ( 2010| report a high 



9 We note that the XMM-Newton study by [Zhang et al.| ( |2006t 
does not classify the system as a merger, underlining again the 
importance of high-resolution X-ray data. 



velocity dispersion of 1300 km s~ for this cluster), we may 
also simply be catching this merger at a later stage than 
MACS J0025. 4-1 222. Several spectacular gravitational arcs 
visible in our UH2.2m images were also independently noted 
by |Hennawi et al.| ( |2008t and ( [Horesh et al]|2010p . An in- 
depth study of the mass distribution in the cluster core is 
underway based on strong lensing constraints derived from 
HST imaging data and Keck spectroscopy. 
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5.3 Redshift evolution of the merger fraction 

Theories of hierarchical structure formation predict an in- 
crease in the fraction of unrelaxed clusters as we approach 
the epoch of formation for massive clusters at z ^ 1. Al- 
though the redshift range of our sample (0.15 < z < 0.7) falls 
short of probing this era we may still expect to see the im- 
print of hierarchical merging in the form of an increase with 
redshift in the fraction of dynamically disturbed clusters, 
as evidenced by higher X-ray/optical separations and by a 
preference for higher values in the assigned morphological 

cod(EE] 

In order to ensure that any observed trends with red- 
shift are not caused by selection bias (not all clusters in our 
original list have Chandra data), all results presented and 
discussed in this section are based on the analysis of a sub- 
set of only 75 clusters, defined by a more aggressive X-ray 
luminosity cut of Lx,rass > 7.5 x 10^^ erg s~^ (and shown 
as a dashed line in Fig.[T]). Above this threshold 95% of the 
clusters in our sample have CXO data, rendering selection 
bias negligible. 

Figure [6] shows the cumulative distribution of the off- 
set of the BCG from the X-ray peak (or centroid) for clus- 
ters above and below the median redshift of this high- 
completeness sample. A two-sided Kolmogorov-Smirnov 
(KS) test finds the probability of the two data sets to be 
drawn from the same parent population to be 44% and 0.3% 
for the BCG/X-ray peak and BCG/X-ray centroid sepa- 
ration, respectively. Fig. [7| shows the KS probabilities for 
the same comparison but as a function of the redshift at 
which we divide the high-completeness sample. Both panels 
of Fig. [T] show a clear increase in the fraction of disturbed 
clusters with redshift, indicating a higher incidence of large 
offsets for the more distant clusters. However, the redshift 
at which this trend becomes significant depends markedly 
on the diagnostic used {z ^ 0.45 for the offset between BCG 
and X-ray peak; z ^ 0.35 for the offset between BCG and 
X-ray centroid). Such a systematic difference does not come 
unexpected though, since the two kinds of offsets probe dif- 
ferent spatial (and temporal) scales. In a post-merger sys- 
tem, a cool core close to either BCG will either have survived 
the collision, or will form as soon as partial relaxation al- 
lows radiative cooling of a stable region of high gas density 
around the BCG. As a result. X-ray/optical alignment will 
be restored quickly on small scales. By contrast, realignment 
between BCG and X-ray centroid after a merger will only 
occur much later, when relaxation has progressed to a global 
scale. 

The same trend is apparent in Fig. [S] which shows the 
redshift dependence of the fraction of clusters with offsets 
in excess of two threshold values, namely the mean, and the 
mean plus one standard deviation, of the X-ray/optical sep- 
aration. Similarly, Fig.[9]probes any variations with redshift 
in the fraction of systems with morphological classifications 
^ 2, 3, and 4. A highly significant increase with redshift is 
observed in the fraction of disturbed clusters, regardless of 
the diagnostic used. 



Note that our inability to reliably identify mergers proceeding 
close to our line of sight does not affect our measurement of the 
evolution of the merger fraction since the distribution of merger 
axes does not correlate with redshift. 



The systems classified by us as the most extreme (i.e., 
the ones featuring the highest X-ray /optical separations or 
a morphological code of 4), however, exhibit little, if any, 
evolution with redshift. This flattening of the redshift dis- 
tribution of the most obvious mergers is likely to be due 
to two effects. The first one is the fact that our classifica- 
tion is based on projected, not three-dimensional separa- 
tions. Hence, systems merging approximately in the plane 
of the sky will be ranked as the most extreme (by all of 
our metrics), while mergers featuring entirely comparable 
physical separations (or morphological complexity) , but pro- 
ceeding along axes closer to our line of sight, may still be 
recognised as disturbed, but will be assigned lower values for 
the (projected) X-ray/optical separation and morphological 
class. This effect causes the redshift distributions of clusters 
of differing degrees of relaxation to scale approximately by 
factors that correspond to different effective inclination an- 
gles of the merger axis. A second effect is clearly at work 
though in Figs, [s] (right panel) and [9] which show a flat- 
ter redshift distribution for the most extreme mergers than 
could be explained by such scaling. Although the still rela- 
tively small size of our sample prevents us from quantifying 
the effect, we suspect that the low X-ray surface brightness 
of the most disturbed systems caused them to be detected 
less reliably near the flux limit of the RASS, resulting in a 
bias against very extended mergers (class 4 clusters) at the 
high-redshift end of our study. The ultimate cause of such a 
bias, namely the combination of poor photon statistics (at 
z > OA the median net photon count in the RASS for our 
sample is 30) with the use of a detection algorithm optimised 
for the detection of point sources, has been identified before. 



Ebeling et al.| ( ^2000b ) , investigating the failure of the Ein- 
stein Medium Sensitivity Survey to detect C1J0152. 7-1357, 
a BHOM z — 0.83, in spite of sufficient photon statis- 
tics in the Einstein raw data, find that complex morphology 
causes the significance of sources near the detection limit to 
be underestimated, resulting in a bias against mergers with 
pronounced substructure. 



6 SUMMARY 

In order to identify the most extreme cluster mergers and 
characterize their evolution with redshift, we have applied 
a set of simple X-ray /optical diagnostics to a statistically 
complete sample of 108 clusters at 0.15 < z < 0.70 and at 
declinations north of —40°, featuring X-ray luminosities of 
Lx > 5 X 10^^ erg s~^ (0.1-2.4 keV), and observed with 
the ACIS imaging spectrometer aboard Chandra. We deter- 
mine the projected (2D) separation between the X-ray peak 
or the X-ray centroid and the BCG, and perform a visual 
morphological classification following the criteria of ^Ebelingl 
et al.|(2007). 



We identify 27 clusters as extreme and active mergers by 
requiring such systems to feature a morphological code of 3 
or 4, and values of at least 71 kpc for the offset between BCG 
and X-ray centroid, or 42 kpc for the separation between 
BCG and X-ray peak. Note that, given the nature of our 
diagnostics, these criteria will favour the selection of mergers 
proceeding along an axis close to the plane of the sky. 

Of those 27, 17 exhibit X-ray/optical characteristics 
suggestive of a binary, head-on merger (BHOM), making 



Cluster classification and evolution 11 



them promising targets for quantitative studies of, e.g., 
galaxy-gas interactions, ICM shocks, or the colUsional prop- 
erties of dark matter. This hst includes MACS J0025.4-1222 
and A 2146, both of which have already proven extremely 
valuable for such work, and also A 1758 which is the target 
of an ongoing multi-wavelength follow-up campaign. Impor- 
tantly, more than half of this list of potential BHOMs are, 
however, yet to be studied in depth. 

The remainder of our initial sample of 27 extreme merg- 
ers are likely to be too complex to be easily modeled, ren- 
dering the physical interpretation of the observables (all of 
them accessible only in projection) difficult at best. Such 
challenges notwithstanding, our list of 10 complex mergers 
still holds great promise for, e.g., statistical studies of the 
properties of gas and galaxies in cluster mergers, as well as 
for investigations of the origin of diffuse, non-thermal ra- 
dio emission which is indeed ubiquitous among the mergers 
selected here ( jCassano et al.|2010| ) . Among this second sub- 
sample to emerge from this study are several famous clus- 
ters, such as A 2744 and the most disturbed cluster presently 
known at z > 0.5, MACS J0717.5+3745. Although many 
systems on this list have been known to be complex merg- 
ers before, about one third are new identifications resulting 
from this work. 

In order to probe the evolution of the merger fraction 
on cosmological timescales, we examine the redshift distribu- 
tion of the clusters in our sample identified as unrelaxed by 
each of our three selection criteria. We find strong evidence 
for such evolution regardless of the merger criterion and re- 
gardless of the statistic used. The increase with redshift of 
the merger fraction starts approximately at z — 0.4; the ob- 
served dependence of this value on the kind of X-ray/optical 
separation used (relative to the peak or the centroid of the 
X-ray emission) is consistent with the difference in relax- 
ation scales (both spatial and temporal) probed by the two 
diagnostics. The fraction of fully relaxed clusters decreases 
dramatically with lookback time even within the limited red- 
shift range of our sample, lending strong support to previous 
claims of a pronounced dearth of cool-core clusters at high 
redshift (fVikhhnin et al]|2007| [Santos et al . 2008; SamueTej 
et al. 2011). Since our analysis is based solely on imaging 



data and thus on quantities observed in projection, mergers 
proceeding at angles close to our line of sight have a high 
probability of being missed. As a result, the true evolution 
of the merger fraction with redshift is almost certainly even 
stronger than reported here. 

Finally, 24 clusters are highlighted in Table 1 as new 
releases from the Massive Cluster Survey, bringing the total 
to 70. We note that, contrary to the previous MACS releases 
( Ebeling et al.|2007[f2010 ), this set of 24 does not represent 
a statistically complete subsample. 



7 OUTLOOK 

Deeper Chandra observations of several clusters in our sam- 
ple with currently poor photon statistics will allow a more 
robust determination of the merger fraction around z — OA 
and thus a more accurate measurement of cluster evolution 
at this dynamically important juncture. In addition, extend- 
ing this sample of X-ray luminous clusters to higher red- 
shift could provide significant constraints closer to the era 



of cluster formation, especially with regard to the evolution 
of cool-core clusters. Finally, we aim to overcome the pri- 
mary limitation of the study discussed here (an unavoidable 
bias against mergers proceeding along an axis close to our 
line of sight) by including dynamical information (ICM gas 
temperatures and radial velocity dispersions of the cluster 
galaxy population) that will allow us to extend the identifi- 
cation of cluster mergers to three dimensions. 
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Figure 1. Luminosity-redshift distribution of all 129 clusters 
meeting our initial selection criteria (see Section|3|. Black symbols 
mark the 108 clusters that have been observed with CXO/ACIS 
and are thus included in this study, red symbols represent the 21 
clusters without CXO data. The X-ray luminosity values corre- 
spond to the flux enclosed in the RASS detect cell. The overall 
sample is volume limited out to z ~ 0.4 and becomes flux lim- 
ited thereafter. Less X-ray luminous lusters are less likely to have 
been observed by CXO. The dashed line marks the limiting X- 
ray luminosity of a nearly complete subsample for which the CXO 
coverage is 95%. 
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Figure 2. VRI image of A 2744 obtained with the UH2.2m telescope; overlaid are logarithmically spaced isodensity contours of the 
adapt ively smoothed X-ray surface brightness as observed with CXO/ACIS-I. The X-ray contours in this and all other X-ray /optical 
overlays (see figure captions for exceptions) are placed at (1.2 x 10~^) x 1.4^ counts s"-*^ in the 0.5-7 keV range. The right panel shows 
a combined exposure, with a total of 74 ks (resulting in over 36,000 counts). The left panel is created by scaling the effective exposure 
time down (1.5 ks exposure time, resulting in 1100 net counts) to match the total counts from the poorest CXO data in our sample (see 
Fig. [3}. The comparison of the two overlays shown here illustrates that even dramatic differences in CXO/ACIS exposure time have 
little impact on our quantitative and qualitative diagnostics (BCG-X-ray peak, and BCG-X-ray-centroid separation, as well as overall 
morphological classification). Note, however, that in this particular case the location of the peak of the X-ray emission shifts significantly 
from the approximate centroid of the southeastern (main) cluster to the cold core almost due south of it (|Owers et al.|2011|. 
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Figure 3. Total number of X-ray photons (after background sub- 
traction) for each cluster within the lowest X-ray contour shown 
in the respective X-ray /optical overlay as a function of redshift. 
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Figure 4. Histogram of astrometric corrections applied to the 
imaging data for our 108 clusters. The median correction is 0.36^^ 
for all optical sources. Optical data from the UH2.2m telescope 
require considerably larger corrections (median ~0.8^^) than data 
from the SDSS and DPOSS2. Three anomalous cluster fields that 
require astrometric corrections in excess of 1.2" are not shown. 
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Figure 5. Histogram of the offsets of the BCG from the X-ray 
peak (top) and X-ray centroid (bottom) for all 108 clusters in our 
sample. The distributions are lognormal centered at 11.5^g° and 
21.2^^g for the peak and centroid respectively. Poisson errors are 
shown, as are the cumulative distributions of offsets. 
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Figure 6. Cumulative distribution of the offsets of the BCG from 
the X-ray peak (top) and X-ray centroid (bottom) for two subsets 
of our high-completeness sample, comprising clusters with red- 
shifts below and above the median redshift of ^=0.387 (solid and 
dotted lines, respectively). While the distribution of offsets be- 
tween BCG and X-ray peak (left) changes little between subsam- 
ples generated by this particular redshift cut, we see a pronounced 
shift to larger offsets of the BCG from the X-ray centroid for the 
more distant clusters (right). A Kolmogorov-Smirnov test yields 
probabilities of 44 (0.3)% for the low- and high-redshift distri- 
butions of the observed BCG-X-ray-peak (BCG-X-ray-centroid) 
values being drawn from the same parent distribution. 



Figure 7. In order to quantify any systematic redshift depen- 
dence of the X-ray/optical offsets, we split our high-completeness 
sample at a given redshift and compare the offset distributions 
of the two subsamples. Shown here are the Kolmogorov-Smirnov 
probabilities for the resulting distributions of BCG/X-ray peak 
offsets (top) and of BCG/X-ray centroid offsets (bottom) as a 
function of the redshift used to divide the sample. We use a run- 
ning median of the Kolmogorov-Smirnov probability to reduce 
statistical variations (although they are still present), and require 
at least 7 clusters in either subsample, which leads to the effective 
exclusion of the shaded redshift range. In each panel the dashed 
line marks the median redshift used to create the subsets com- 
pared in Fig. [6] 
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Figure 9. The fraction of clusters with morphological classifica- 
tions ^ 2 (black), ^ 3 (blue), and = 4 (red). Again the fraction 
of disturbed clusters (all systems except the ones classified as 
fully relaxed) clearly increases with redshift (significant at the 
5.0(7 confidence level), consistent with the trends apparent from 
Figs. [T] and [S] See text for further discussion. 



Figure 8. Top: The fraction of clusters with offsets of the BCG 
from the X-ray peak in excess of 12kpc (the mean of the dis- 
tribution of offsets; black), and 42kpc (mean plus one standard 
deviation; red) as a function of redshift. The evolution visible in 
the black histogram is significant at the 2.7a confidence level. 
Bottom: Same as left panel, but for the offsets of the BCG from 
the X-ray centroid, and using thresholds of 21 (mean) and 71 kpc 
(mean plus Icr). The evolution observed for the larger sample is 
significant at the 4.5cr confidence level. For both panels the bin 
boundaries are chosen such that each bin contains 15 clusters; 
errors are computed following |De Propris et aL\ ( |2004| ) . The evo- 
lutionary trend apparent in either panel and the approximate red- 
shift at which the merger fraction increases significantly are fully 
consistent with the results shown in Fig. [T] based on unbinned 
data. 
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" BHOM or complex merger; see Figs. |10|11|12| 

^ Cluster detected by the Massive Cluster Survey (MACS). 

Previous releases of MACS clusters are presented and dis- 



cussed in Ebeling et al. (2007 2010) 



Note. -All luminosities as measured from ROSAT All-Sky 
Survey data within the nominal detect cell. Morphological 
classification codes are explained in Section [2] 



Referen ces for Table 1 (1 ) This work , (2) [Allen et al 



(7) Ebehng et al. (2000a|), (8) 



Ebeling et al. 



(1992), (3) Bohringer et al. (2004), (4) Crawford et al. 



( T995t, (5)|Crawford et al.l ( |1999|), (6)|Ebeling et al.|(|1998| ). 



Ebeling et aLn2007), (9) 



2010| ), (10) [Elvis et al.| ( |1992t , (11) [Gioia 

Luppino|([1994t , ( 12) [Le Borgne et al] |l99l| , (13) [Oegerle 
et al.[([1991[), (14) [Romer et al.[ (|1994|), (15) Sarazin, Rood 



k Struble' f 19821 , (16) |Stocke et al.[ " p99Tt, (17) [Struble k 
Rood ( 1999J , (18) [White| ( |2000[ ), (19) [Wright, Abies Allen 
( T983| . 
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Figure 10. Logarithmically spaced isodensity contours of the adapt ively smoothed X-ray surface brightness as observed with Chandra 
overlaid on optical images of the clusters identified as complex mergers (see Section [5. 2. 1| . The units on the X-ray contours are (1.2 x 
10-6) X 1.4^ photons s-i in the 0.5-7 keV range. Shown are (in R.A. order) A 2744, A 520, MACS J0717. 5+3745, A 665, A 1300, 
MACS JII49. 5+2223, A 1682, MACS J2129.4-0741, RX J2228. 5+2036, MACS J2243. 3-0935, and A 2631. High resolution versions of 
each overlay is available at http://ifa.hawaii.edu/~amann/MNRAS2011highres/ 
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Figure 11. Same as Fig. [Tojfor the primary BHOM candidates discussed in Section [5. 2. 2| Shown are (in R.A. order) MACS J0025.4- 
1222, MACS JOI40. 0-0555, MACS J04 1 7.5-1 154, A 520, MACS J0553. 4-3342, MACS J1006. 9+3200, A 1758, A 1914, A2146, and 
MACS J1731. 6+2252. 
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Figure 12. Same as Fig. [To] for the secondary BHOM candidates. Shown are (in R.A. order) A 2813, MACS J0358. 8-2955, 
MACS JO404. 6+1109, MACS JO4I6. 1-2403, ZwCl 0847.2+3617, MACS J1115. 2+5320. 
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Right Ascension (2000) 






56" 54' 52" 



Rigiit Ascension (2000) 






Right Ascension (2000) 



46" or 59" 4 




Right Ascension (2000) 




Right Ascension (2000) 




Right Ascension (2000) 



Figure Al. (Continued over next 6 pages). Same as Fig. |10| for the all clusters in our sample not shown in Fig. |10| or |12| Shown 
are (in R.A. order) MACS 001 1.7-1523, CL 0016+1609, MACS J0035.4-2015, ZwCl OIO4+OO48, MACS JOlll. 5+0855, A 209, RBS436, 
MACS J0159. 0-3412, MACSJOI59.8-O849, MACSJO242.6-2132, MACS J0257. 1-2325, and A 402. High resolution versions of each 
overlay is available at http: / / ifa.hawaii.edu / ^amann/MN RAS2Ql Ihighres / 1 
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Figure Al. continued Shown are (in R.A. order) A 3088, MACS J0308. 9+2645, MACS J0326. 8-0043, MACS J0329. 6-021 1, 
MACS JO429. 6-0253, RX J0437. 1+OO43, RX J0439. 0+0715, MACS JO451. 9+0006, MS 0451.6-0305, MACS J0455. 2+ 0657, and 
MACS J0520. 7-1328, and A 907. 
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Figure Al. continued Shown are (in R.A. order) MACS J0547.0-3904, MACS J0647.7+ 7015, A 586, ZwCl 0735.7+7421, 
MACS JO744. 8+3927, A 661, A 697, MACS J0911. 2+1746, MACS J0913. 7+4056, A 773, RBS 0797, and ZwCl 0947.2+1723. 
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Figure Al. continued Shown are (in R.A. order) ZwCl 1021.0+0426, A 3444, MACS Jl 105.7-1014, MACS J1108. 8+0906, 
MACS J1115. 8+0129, MACS J1206. 2-0847, MACS J1226. 8+2153, MACS J131 1.0-0310, A 1689, A 1722, A 1763, and RX J1347. 5-1145 
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Figure Al. continued Shown are (in R.A. order) MACS J1359. 1-1929, A 1835, MACS JI4II. 3+5212, MACS JI423.8+2404, 
MACS JI427. 2+4407, MACS JI427. 6-2521, MACS JI447. 4+0827, A 1995, ZwCl 1454.8+2233, AS780, ZwCl I459.4+424O, and 
RBSI46O 
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Figure Al. continued Shown are (in R.A. order) RX J1532. 9+3021, MACS J1621. 3+3810, A 2204, A 2219, RXC J1720. 1+2637, 
MACS J1720. 2+3536, A 2261, MACS J1931. 8-2634, RXC J2014.8-2430, MACS J2046. 0-3430, MACS J2049. 9-3217, and 
RXJ2129. 6+0005 
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